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Abstract 

The  purpose  of  this  report  is  to  communicate  recent  results  on  the 
existence  and  uniqueness  of  solutions  to  the  Vlasov  equation.   Both  the 
Vlasov-Poisson  system  and  the  Vlasov-Maxwell  systems  are  considered. 
Chapter  1  is  devoted  to  the  Vlasov-Poisson  system.   A  theorem  is  stated 
which  gives  conditions  for  the  solvability  of  the  system  for  general 
once  dif ferentiable  initial  data  having  compact  support.   As  an  example 
of  the  application  of  the  general  theorem  the  case  where  the  system  has 
spherical  symmetry  is  considered.   A  global-in-time  existence  and 
uniqueness  theorem  is  proved  for  this  case.   Chapter  2  deals  with  the 
Vlasov-Maxwell  system.   Here  a  local-in-time  existence  and  uniqueness 
theorem  is  proved  for  a  general  class  of  initial  data  having  compact 
support. 
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Chapter  1:   The  Viasov-Poisson  system 

I.   Introduction 

The  Viasov-Poisson  system  of  equations  describes  the  motion  of  a 
collection  of  particles  in  the  presence  of  its  own  electrostatic  (or 
gravitational)  field.   The  system  under  consideration  is 

a)  |4  +  Wxf  +V^^^   =  0 

(1.1)  f(x,v,0)  =  g(x,v) 

b)  At  =  -4tt  2  /  (f-b)dv 

(x,v)  -  a  point  in  6-D  phase  space 

f(x,v,t)  -  phase  space  distribution  function 

b(x,v)  -  fixed  background  distribution 
We  assume  g,b  Ci(R^),  i.e.,  continuously  dif f erentiable  functions  with 
bounded  support.   The  support  of  a  function  is  the  region  of  phase 
space  outside  of  which  the  function  is  zero.   A  function  with  bounded 
support  goes  to  zero  outside  a  bounded  region  of  phase  space.   In  this 
paper  a  general  theorem  is  stated  which  gives  conditions  for  the 
solvability  of  system  (1.1).   The  proof  of  the  theorem  is  given  in 
[14].   The  bibliography  at  the  end  of  this  chapter  contains  a  set  of 
references  for  the  classical  existence  theory  of  (1.1).   An  application 
of  the  general  theorem  is  made  to  the  case  where  the  initial  data  and 
hence  the  solution  have  spherical  symmetry.   The  result  is  a 
global-in-time  existence  and  uniqueness  theorem  for  the  system  with 
spherical  symmetry. 
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II.   General  theorem 

The  following  notation  is  used 
R^  -  six  dimensional  phase  space 
Cq(R^)  -  the  space  of  continuous  functions  with  continuous 

bounded  derivatives  on  Rr   and  which  have  bounded  support  in  R/-. 
C  (R^x[0,T])  -  The  space  of  continuous  functions  with  continuous 

bounded  derivatives  on  R/-x[0,T  ]  ,T-time. 

The  term  "solution"  refers  to  a  function  that  is  of  class  C  , 
satisfies  (1.1),  and  which  is  integrable  over  phase  space  in  such  a  way 
that  appropriate  conservation  laws  (like  mass,  momentum,  energy)  can 
hold. 

Theorem:   Let  g,b  be  functions  in  C^(R6).   The  solution  to  (l.la.b) 
exists  and  is  unique  as  an  element  of  C  (R^x[0,T])  if  and  only  if  the 
system  admits  an  a  priori  bound  on  the  support  of  C*  solutions  for  t  in 
the  interval  [0,T] . 

Proof:  see  [14]. 

This  theorem  gives  a  criterion  to  determine  whether  the  solution 
to  (l.la,b)  exists  and  is  unique  for  a  particular  initial 
configuration.   To  use  the  theorem  let  us  make  the  assumption  that  a  C 
solution  to  the  problem  exists  on  the  interval  of  time  [0,T].   The 
theorem  tells  us  that  a  C1  solution  to  the  problem  necessarily  has 
bounded  support  in  phase  space.   We  then  try  to  compute  a  priori  what 
the  bound  on  the  support  of  such  a  solution  must  be  as  a  function  of 
time.   If  it  is  possible  to  compute  such  an  a  priori  bound  for  a 
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solution,  then  the  theorem  teiis  us  that  the  C1  solution  to  the  system 
in  fact  exists,  is  unique,  and  of  course  has  the  bound  on  support 
already  computed.   The  important  bound  to  compute  is  on  velocity. 
For  initial  data  with  bounded  support  the  statement  that  the 
support  of  the  solution  remains  bounded  is  equivalent  to  the  statement 
that  the  electric  field  remains  bounded.   Therefore,  an  alternative  way 
of  stating  the  general  theorem  is  as  follows:  a  C   distribution 
function  that  is  a  solution  to  (1.1)  cannot  give  rise  to  an  unbounded 
electric  field  in  a  finite  time.   To  prove  existence  and  uniqueness  of 
the  solution  to  (1.1)  of  class  CVR^xfO.T])  it  is  sufficient  to  show  "a 
priori"  that  the  electric  field  remains  bounded  for  t  e  [0,T]. 


III.   The  spherically  symmetric  Vlasov-Poisson  system 

An  existence  and  uniqueness  theorem  for  the  spherically  symmetric 
case  was  first  proved  by  Batt  in  [2],  The  proof  in  [2]  is  valid  for 
the  gravitational  case  alone  where  the  spherically  symmetric  field 
always  points  inward  along  a  radius.   In  [13]  a  different  proof  is 
given  which  makes  no  use  of  an  inward  directed  field  and  is  thus  valid 
for  both  the  gravitational  and  electrostatic  case.   The  proof  in  [13] 
builds,  however,  on  some  of  the  development  in  [2],   In  this  paper  we 
give  an  independent  development  of  the  result  in  [13]. 
Let 


r  =  (xf  +  *l  +  x§)1/2 
u  =  (v2  +  v\  +  v§)1/2 
C  =  xjv^  +  X£V2  +  X3V3 
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Let  the  functions  g,b  in  (l.la,b)  be  of  ciass  Cq(R^)  and  have  the  form 


(3.1)      g(x,v)  =  G(r,u,C)  ,  b(x,v)  =H(r,u,?) 

Initial  and  background  distributions  of  this  form  are  termed 
spherically  symmetric,  and  solutions  to  (1.1)  with  data  of  this  form 
reduce  to  functions  in  the  variables  r,u,5,  and  t  and  are  termed 
spherically  symmetric.   Spherically  symmetric  solutions  are 
characterized  by  the  fact  that  the  field  remains  a  function  only  of 
radius  and  time.   Also  if  (x(t),  v(t))  is  a  trajectory  of  (1.1)  then 

x(t)  x  v(t)  =  c  . 

These  facts  will  become  evident  as  we  proceed. 
Let 

n  =  (xlVl  +  x2v2  +  x3v3)/ru  =  C/ru  . 
For  points  (x,v)  e  R^  such  that  |n  |  t    1  (i.e.  ,  the  velocity  vector 
does  not  point  along  the  radius)  spherically  symmetric  solutions  to 
the  Vlasov  equation  can  be  written  in  terms  of  the  variables 


r  =  (x^  +  x|  +  x^)1/2 


u  =  {v\   +  v\  +   v§)1/2 

a  =  cos   ((xjv.  +  x2v2  +  x-jv3)/ru) 


Here  a  is  the  angle  between  the  position  vector  and  the  velocity  vector 
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and  takes  on  values  between  0  and  it.   In  terms  of  these  variables  the 
Vlasov-Poisson  system  reduces  to  the  form 


N   3F  ,         9F         n,         ,  3F 

a)  - —  +  cos  a  u  - —  -  cos  a  0(t,r)  - — 

3 1  3  r  3u 


3  F 
+  (0/u  -  u/r)  sina  —  =  0 

3a 


(3.2)     F(0)  =  G 


b)  0(t,r)  =  -  ^-  j'   h(t,F)  ?2  dr 
r   0 

h(t,r)  =  2tt  /   /   (F-H)  sina  da  u2  du 
0   0 


where 

H(r,u,a)  is  a  fixed  background  distribution 
h(t,r)  the  charge  density, 
-0  the  radial  electric  field 
The  characteristic  equations  for  (3.2)  are 

r  =  u  cosa 
(3.3)     u  =  -  0(t,r)  cosa 

a  =  (0/u  -  u/r)  sina 

An  integral  of  (3.3)  is 
(3. A)     r(t)u(t)  sin(a(t))  =  r(0)u(0)  sin(a(0))  =  c  . 

Thus  if  a  trajectory  of  (3.2)  (solution  to  (3.3))  has  its  initial  point 
in  the  region  n  *    1  then  the  trajectory  remains  in  the  region  n  *  1. 
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That  is  n  M  implies  r(0)u(0)  sin(ct(0))  =  c  *  0.  From  (3.4)  it  follows 
that  r(t)u(t)  sin  (a(t))  *  0  for  all  t  and  hence  the  trajectory  remains 
in  the  region  n  *  1  for  ail  t.  In  order  to  prove  the  result  of  giobai 
existence  and  uniqueness  for  the  spherically  symmetric  Viasov-Poisson 
system  we  analyze  trajectories  in  the  region  n  *  1  and  compute  bounds 
on  the  trajectories  in  this  region.   By  continuity  results  for  the 
region  n  t    1  are  extended  to  the  lines  n  =  1.  By  this  means  an  a  priori 
bound  on  velocities  in  the  system  is  computed. 


IV.   Bounds  on  the  field 

Let  the  initial  data  G  to  (3.2a,b)  have  the  bound 

0  <  G  <  D 
and  G  and  H  have  bounds 

8tt2/   /   J   Gsinada  u2du  r2dr  <.    M/2  , 

0     0      0 

oo    as    fl 

8tt2  /   /    /    H  sina    da    u2du    r2dr  <    M/2 
o   o   o 

Hence    the   solution   to   (3.2)   has   bounds 


oo    oo    7f 

0  <    F  <    D      ,8tt2///f   sina    da   u2du    r2dr  <    M/2 
o   o   o 
and   also 

81,2  I      /      /       |F-H|sinada   u2du    r2dr  <    M 
0     0      0 


Let   us   assume   a   constant  K  such   that 
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|h|    <    2tt    /     |  F  —  H  |    sina    da    u2du  <    S 

The   bounds   D   and  M   on  F   are   known   a   priori    from   the   initial   data.      The 
bound  A"  gets    computed   as    part   of    the   proof. 
Bounds    for   the   field  0    are 


|0  |  <    %jC    |h|    F2   d?<    *X   r3/3 
r2  J  o  r2 


<  Ar,    A  =   4/3  tt   A" 

I©  |   <     l/r2    (4ir    Jr    |h|    F2   dr) 

r  J  0    •     i 

<  1/2    (4ir    /       |h|?2   dr)  <    M,    2 

r  J  0  /r 


A  continuous  function  bounding  0  is 

g(r)  =  Ar  0  <  r  <  (M/A)1/3  =  I 

g(r)  =  M/r2  r  >  £ 

The   continuous   anti-derivative   of   g    is 

G(r)    =  1/2   A   r2-3/2   A1/3  M2/3  r  <    I 

G(r)    =   -M/r  r   >    I 

The   bound   on  0    is 

I©  |   <    G'(r) 


V.  Existence  and  uniqueness  theorem  for  the  spherically  symmetric  case 

Theorem:  Let  the  initial  and  background  distributions  g  and  b  to 
(l.la.b)  be  of  class  Cq(R6)  and  be  of  the  form  (3.1)  (i.e.,  spherically 
symmetric).   For  such  spherically  symmetric  data  the  solution  to 
(l.la.b)  of  class  C  (R^  x  [0,T])  exists  and  is  unique  for  ail  T. 

Proof: 

We  state  without  proof  that  for  initial  and  background  data  of 
form  (3.1)  the  solutions  to  (l.la,b)  satisfy  a  system  of  the  form 
(3.2a,b)  in  the  variables  r,u,a,  and  t  in  the  region  n  *  1.  As  has  been 
shown  trajectories  that  originate  in  the  region  n  *  1  remain  in  this 
region.   An  a  priori  bound  on  the  support  of  the  solution  to  (l.la,b) 
is  obtained  by  analyzing  the  trajectories  (solutions  to  (3.3))  in  the 
region  n  t    1.  The  results  are  extended  by  continuity  to  the  lines  n  = 
1. 

Let  us  assume  the  existence  of  a  solution  F  to  (3.2a,b)  on  the 
interval  of  time  [0,T]  which  has  bounded  support  in  the  variables  r,u. 
We  wish  to  compute  a  priori  what  this  bound  must  be.   Let 

Q(t)  =  supp  F(t)  -  support  of  F  in  r,u,ct  space  for  each  t. 

p  =  (r,u,a)  -  a  point  in  r,u,a  space. 

E(p)  =  u  -  the  projection  of  p  on  the  u  axis. 

fl(T)  =  Ufi(t),  t  e  [0,T] 

P(T)  =  su2  E(p)  -  a  bound  on  the  support 
pen(T) 

in  velocity  space  up  to  time  T. 
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Let  the  constant  A(T)  satisfy 

|h|    <    2tt    /     |  F  —   H  |    sina    da    u2du  <    A(T)    for    t  e     [0,T] 

and   set 

A  =  A(T)    =   4/3ir    ff(T)    . 

We   are  assuming 


8tt2  /      /      /       |F  -  Hi    sina    da   u2du    r2dr  <    M 
0      0  J  0 


Let 

G(r)  =  1/2  Ar2  -  3/2  A1/3  M2/3     r  <  I    =    (M/A)1/3 
G(r)  =  -M/r  r  >  I 

The  bound  on  the  field  0  for  t  e  [0,T]  is 
|0|  <  G'(r) 

For  any  two  points  r^,r2  tne  function  G(r)  satisfies 
(5.1)       |G(ri)  -  G(r2)|  <  3/2  M2/3  A1/3 

Let  p(t)  =  (r(t),  u(t),  a(t))  be  a  trajectory  for  t  e  [0,T]  such 
that  p(0)  is  in  the  region  t\   t    1.  The  proof  now  proceeds  through  a 
sequence  of  steps  to  get  the  estimate 
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u2(t)  <    u2(0)    +   6M2/3   A1/3 

.)      If    r)    0   or    r  <    0  on    [t1,t2]C    [0,T] 

then 

|u2(t2)'-  u2(t1)|   <    2    |G(r(t2))   -  GCrCt^) 
as    follows: 


A.  (1/2  u2)    =  u   u  =    -  0    (t,r)   u    cosa 
dt 


=      -0    r 


Hence 


|u2(t2)    -  u2(t1)|   <    2  /         |0r|dt 


tl 


However 


if      V  >    0  |0r|   <    4-  G(r(t)) 

at 


if      r  <    0  |0r|   <    -A  G(r(t)) 

dt 


Thus 
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.2/a.x  _  .2 


uz(t2)  -  uz(t1)|  s  2  |G(r(t2))  -  G(r(tl)) 


b)   At  a  local  maximum  of  r(t) 

u2(t)  <  M2/3  A1/3 

as  follows: 

At  a  local  maximum 

r  =  u  cosa  =  0  +  a  =  tt  /2 

sin(a)  =  1 
r  =  u  cosa  -  u  sina  a  =  -  u  a  <  0 
->■  a  >  0 

Hence  from 

a  =  [0 /u  -  u/r]  sina  we  have 

0/u  -  u/r  >  0  +   u2  <;  0r 
Thus  0  >    0  i.e.  the  field  -  0  <  0.   From  the  bounds  on  |0 


u2  <:  Ar2      r  <  I 
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uz  <  M/r      r  >   I 


The  maximum  occurs  at  r  =  I   and 


u2(t)  <  M2/3  A1/3 

c)  At  a  local  minimum  of  r(t) 

u2(t)  <  u2(0)  +  3M2/3  A1/3 

as    follows:      Let   t  >    0  be   a   local  minimum  and   suppose   r(t)    =  0.    From 
(3.3)  we  have   that 

r  =  u    coscc    -  u    sin(a)a 
=  -  <J>    cos2a    -  u    sina    [<}> /u   -  u/r]    sina 


2 

9  9  U  9 

4>    cos^a    -  <f>    sin^a   +  —  sin  a 

r 


2 

-i-  u        -2 

+  —  sin  a 


At    a  minimum  a    =  it  /2   and   assuming    r  =   0  then 


u2   =  <br 


Thus    from  the   bounds   on  § 
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u2  <  Ar2     r  <■    I 


uz   <  M/r      r  >   I 


and  ietting  r  =  I    it  follows  that 

u2(t)  <  M2/3  A1/3 

The  estimate  therefore  holds  if  r(t)  =  0. 

Assume  that  r(t)  t  0.  Then  there  exists  tg  <  t  such  that  either 
tQ  =  0  or  tg  is  a  local  maximum  and  r  <  0  on  [tg,t].  If  tQ  =  0  then 
from  (5.1)  and  part  a) 

u2(t)  <  u2(0)  +  3M2/3  A1/3 

If  tg  is  a  maximum  we  distinguish  two  cases: 
i)   r(t),r(tQ)  k    i    or    r(t),r(tQ)  >  I    in  which 

|G(r(t))  -  G(r(tQ))|  <  M2/3  Al/3 
and  from  parts  a,b  and  (5.1) 
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.2(t)  <    .*(*„>  +  *»3  Al/3  «    3M2/3  Al/3 
li)      r(0  <   t,    r(t0)  >    i   in  »hich  from  (5.1)  and   parts  a,b 
At)«    u2(t0)+3M^Al/3<    rf"*"1 

We    further    refine    this    estimate. 
(5.2)  u2(t)  <  u2(t0)   +  2|G(r(t))   -  G(r(t0))| 

<u2(t0)+3M2/3A^-^L_-Ar(t)2 

Since    r(tQ),    r(t)   are  extrema   then 

a(tQ)    =a(t)    =tt/2   and   from   (3.4) 

u(tQ)    r(t0)    =  u(t)r(t) 

r(tQ)  =  r(t)u(t)/u(t0) 
Substituting  in  (5.2)  and  minimizing  the  negative  term  over  r(t)  we  get 


u(tn)2 


u2(t)<u2(to)+3M2/3  .1/3.^/3.1/3^  1/3 
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But  it  is  shown  above  that 


u2(t)  <  AM2/3  A1/3,  thus 


(5.3)     u2(t)  <  u2(tn)  +  JA2'^   A1'3  -  3/4l/3  M2'3  A1/3   (  "^  )U3 

VM2/3A1/3' 


It  is  easiiy  seen  that  for  u  (tg)  in  the  interval  [0,M  '  A   '     ]    (part  b) 

o 

the  right  side  of  (5.3)  is  maximized  at  u  (tg)  =  0.  Thus 


u2(t)  <  3M2/3  A1/3 

This  estimate  is  obtained  in  [2]  through  a  different  derivation, 
d)   For  any  t  >  0 

u2(t)  <  u2(0)  +  6M2/3  A1/3 

as  follows:  If  r(t)  =  0  then  as  before 

u2(t)  <  M2/3  A1/3 

and  the  estimate  holds.  Assume  therefore  that  r(t)  ±  0.  If  t  >  0  there 
exists  tg  <  t  such  that  either  tg  =  0  or  tg  is  a  local  extremum  of  r(t) 
and  r  does  not  change  sign  on  [tg,t].   Thus  from  (5.1)  and  parts  a-c 

5.4)      u2(t)  <  u2(tg)  +  2|G(r(t))  -  G(r(tg))| 
<  u2(0)  +  6M2/3  A1/3 
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Estimate  (5.4)  is  now  used  to  get  an  a  priori  bound  on  the  support 
of  F  in  velocity  space.   In  terms  of  the  bound  on  velocity 

P  =  P(T)   for   t  e  [0,T] 
the  charge  density  h  has  the  bound 


00       "FT  o 

|h|    <    2tt    /      (       IF   -  Hi    sina    da   uzdu 
0  ;  0 

<    !Z.  D   P3(T)    +  B   <    !?L  (D  +   B)P3(T)   <    ^-  CP3(T) 


(assuming   P(T)   >    1).      Here  D   =  sup|F|    , 

oo    it 

2tt    /   j     |H  |    sina    da   u2du  <    B    ,    C  =   D  +  B 
o  o 

Therefore   let 


X(T)    =  ^L  C  P3(T) 


and   set 


A   =  A(T)    =   4/3  tt    A(T)    =   (-^-)  2   C   P3(T) 


=  (^-)2c  P3 


From  (5.4)  it  then  follows  that 


(5.5)      u2(t)  <  u2  +  6M2/3  {^-)2^   C1'3  P 
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Thls  Inequality  Is  satisfied  for  any  solution  to  (3.3).   In  particular 
it  Is  satisfied  for  p(0)  e  fl(0),  the  support  of  the  Initial  data  G. 
Since 

P  -  sup  u(t),  t  e  [0,T] 
p(0)en(0) 

It  follows  from  (5.5)  that  an  inequality  for  P  is 

P2<  P2  +  6M2'3  (^-)2/3  C1'3  P 

PQ  =  supE(p) 
pefi(O) 

Thus  an  inequality  for  P  is  of  the  form 


P2  <  Pq  +  B  M2/3C1/3P 


where  M  and  C  are  constants  known  a  priori  from  the  initial  and 
background  distributions.   A  bound  for  P  is 


P  =  P(T)  <  PQ  +  B2M2/3C1/3 


The  bound  is  independent  of  T  and  is  therefore  uniform  for  all  T.  The 
uniform  bound  on  support  in  velocity  space  leads  to  a  bound  on  support 
in  position  space.   From  the  general  theorem  we  thus  get  the 
global-in-time  existence  and  uniqueness  of  the  classical  solution  to 
(l.la.b)  with  data  of  the  form  (3.1). 
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Chapter  2:   The  Vlasov-Maxwell  system 
I.   Introduction 

The  purpose  of  the  present  paper  is  to  give  a  local  exis- 
tence and  uniqueness  theorem  for  the  Vlasov-Maxwell  system  cf 
equations.   The  Vlasov-Maxwell  system  is  one  of  the  fundamental 
equations  underlying  the  kinetic  theory  of  plasma.   It  is  a 
non-linear  integro-dif f erential  system  which  describes  the 
motion  of  a  collection  of  charged  particles  in  the  presence  of 
its  own  self  generated  electro-magnetic  field.   Under  discussion 
is  the  following  initial  value  problem 


a)    |t+W  f-  (E  +  vx  b)«V  f  =  0,    f  (0)  =  g 
3t       x  v 


(1.1) 


b)    ||  -  V  x  B  =  4tt 


vf  dv  ,  ~  +  V  x  e  =  0 

o  t 


where 


E(0)  =  E  .  B(0)  =  B^ 
o  o 


c)    V'E   =  -  4tt 
o 


gdv  ,  V • B   =0 
o 


E  =  (E1,E2,E3)  ,  B  =  (B1#B2,B  ) 


x  =  (x1,x2,x3)  ,  v  =  (v1'v2'v3) 


The  function,  f,  gives  the  number  density  of  particles  in  phase 
space.  The  functions  E,  B  are  the  electric  and  magnetic  fields 
produced  by  the  distribution,  f.  The  system  is  internally  con- 
sistent in  that  the  distribution,  f,  evolves  in  time  under  the 
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influence  of  the  fields  E,  B,  which  are  in  turn  being  generated 
by  f.   Given  initial  condition  (1.1c)  it  is  easy  to  see  that 
the  solution  (1.1)  satisfies  the  additional  equations 

V«E  =  -  4tt    fdv  ,  V«B  =  0  for  t  >  0  . 


In  the  present  context  f  represents  a  single  species  of  charged 
particles,  say  electrons.   Normally,  a  plasma  is  comprised  of 
multiple  species  each  satisfying  a  transport  equation  of  the 
type  (1.1a)  and  each  contributing  to  the  charge  and  current 
terms  on  the  right  side  of  equation  (l.lb,c).   The  analysis 
given  in  this  paper  can  be  modified  to  deal  with  such  an  expanded 
multiple  species  system  of  equations.   For  the  sake  of  simplicity 
we  deal,  however,  with  the  single  species  system  given  by  (l.la,b,c) 

The  existence  and  uniqueness  theorem  for  (1.1)  is  produced 
by  generalizing  a  theorem  of  Kato  [  3,  theorem  II,  p.  195]  .   In  [  3] 
Kato  proves  a  local  existence  and  uniqueness  theorem  for  quasi- 
linear  symmetric  hyperbolic  systems  of  equations.   In  the  proof 
the  author  actually  generalizes  somewhat  the  type  of  systems  being 
considered.   The  coefficients  in  the  equations  are  not  merely 
point  wise  operators  on  functions  of  x  and  t  but  are  instead 
operators  which  are  point  wise  in  t  but  which  are  non-local  in 

x  €  R   mapping  functions  of  R   into  functions  on  R  .   In  the  pre- 
n  n  n  r 

sent  paper  we  further  generalize  the  result  in  [  3]  so  that  the  co- 
efficients are  non-local  in  both  x  and  t  mapping  functions  of 
R  x[  0,T]  into  functions  on  R   x[  0,T]  .   The  starting  point  of 
this  analysis  is  the  linear  theorem  [3,  Theorem  I,  p.  189]  for 
linear  symmetric  hyperbolic  systems.   This  theorem  is  used  as 
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in  [  3]  to  prove    local  existence  and  uniqueness  for  a  general  sys- 
tem with  operator  coefficients  of  the  type  mentioned  above. 
Solutions  are  obtained  as  continuous  mappings  of  t  into  classes  of 
Sobolev  spaces  on  R  .   The  proof  given  is  in  fact  a  modification 
of  that  in  [  3]  .   Having  generalized  Kato's  work  we  can  then  give 
a  proof  for  local  existence  and  uniqueness  for  the  Vlasov-Maxwell 
system.   The  proof  we  give  is  for  an  initial  distribution  with 
compact  support.   For  this  class  of  data  the  Vlasov-Maxwell  sys- 
tem can  be  put  into  a  form  so  that  the  general  theorem  applies. 

As  far  as  we  know  the  existence  and  uniqueness  theorem  pre- 
sented in  this  paper  is  the  first  such  result  obtained  for  the 
fully  non-linear  three-dimensional  Vlasov-Maxwell  system.   Pre- 
viously some  results  have  been  obtained  for  the  system  in  lower 
dimensions  and  with  various  additional  modifications  or  restric- 
tions, [1],  [  2]  ,  [  4]  . 
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II.   Function  spaces  and  notation 

We  are  adopting  as  much  as  possible  the  notation  of  [3],  Sore  chanaes 
in  notation,  however,  are  made. 


Rj^  -  m  dimensional  Euclidean  space 
P  -  A  real  or  complex  Hilbert  space  which  in  the  present  context 
will  be  identified  with  ^  for  some  value  of  m. 


m 

For  x  6  IV  ,  |x|  =  y^Xi2)1/2 


u(Rm)  -  a  real  P  valued  functions  defined  on  ^ 

L^R^P)  -  the  linear  space  of  P  valued  functions  on  R^   such  that 


u|lo,m=  (  1  |u|2dx)l/2  < 


Hs(Rm,p)  -  the  space  of  P  value  functions  on  ^  for  which  the 

distributions  derivatives  of  order  <_  s  are  of  class  L  (R_,P). 


The  norm  is 


UHs,m  =  U,     Ca||tf<u||2/m)l/2 


ai  <s 
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Ca  -  positive  constants 


Lj(R_,p)   -  the  set  of  P  valued   functions  on  R^  which  are 
uniformly  locally  L2  in  the  sense  that 


HuHo£,m  =  sup     [I  |u(y)|2dy]1/2    < 

xe^  |  x-y|  <  1 


^(R^P)   -  the  set  of  P  valued  functions  on  P^  whose  distribution 

derivatives  of  order  _<_  s  are  of  class  L^CR^p).     The  norm  is 


llull  s£,m  =  sup<s||EPu||oJl/m  . 


Let  B(P)  be  the  3et  of  bounded  linear  operators  on  P.     We  are  only  interested 
in  the  case  where  P  is  finite  dimensional  in  which  case  elements  of  B(P)  are 
pxp  matrices   (p-dimension  of  P).      For  a  =  [a^    .]  e  B(P)  a  suitable 
operator  norm  is 


1*1 -&     |ai-i 

-1-'  J 


|2ll/2 


1^(1^, B(P))  -  the  set  of  functions  of  T^  valued  in  B(P)  for  which 
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a||oJ,  _  =  suo  ;a(y)|  2dy'  l/2   < 


H|(R^B(P))   -  the  set  of  functions  on  F^  valued   in  B(P)  whose 
derivatives  of  order  _<_  s  are  elements  of  Lf  (R^BCP) ) . 


^  II  -  f V       r    II  rPa  II2       I1/2 

a  Ils4,m  "  U<=     cH  NoA,mJ 


i  ot  I  <s 


C  -  postive  constants 


L^CCT.'H^R^P)]  -  the  set  of  P  valued  integrable  functions  u(t,«  ) 
on  [O.Tlxl^  such  that  for  each  t  u(t,' )  €  H^R^P). 

CntO,T;Hs(Rm,P)]  -  the  set  of  P  valued  function,  u(t,0,  on 

[OrTlxl^  such  that  the  functions  Dau(t,*  ),  |a|_<  m,  are 

s 
continuous  mappings  of  t  into  H  (R  ,P). 

Similar  definitions  apply  for  mappings  of  t  into  H|(R^,P), 
H^CRj^BfP)).  As  normed  the  spaces  H^R^P)  are  Hilbert  spaces 
^(Rm/P)/  Hf^'8^))  are  Banach  spaces. 


-26- 


III.   Linear  theorem 

In  this  section  we  state  the  linear  theorem  on  which  the 
proofs  for  the  non-linear  results  depend,  this  theorem  is  proved 
in  [3] .   The  proof  is  rather  complex  and  is  an  application  of  the 
theory  of  semigroups   to  evolution  equations.   We  will  not  give 
the  proof  here  but  will  simply  state  the  result.   For  simplicity 
a  somewhat  less  general  theorem  is  stated  but  one  which  is  ade- 
quate for  our  purposes. 

The  following  Cauchy  problem  is  considered 


a  m  „ 

(3.1)  |H  +      i      a    (t,x)    |£-  =   f (t,x) 

j=l      ■>  j 


0<t<T,x6R      ,    u(0,x)    =   u    (x) 
—        —  m  o 


the  functions  u  ,  f  are  valued  in  the  Hilbert  space  P  which  in  the 

present  context  is  taken  to  be  R   for  some  n  ,  the  coefficients  a  . 

c  n             '                                                j 

are  in  B (P) . 

Theorem  3.1.  Let  s  be  an  integer  such  that  s  >  m/2  +  1  ,  and 

■ 

let  1  <  s   <  s  .  Assume  the  following  conditions: 


(3.2)      a..  6  C[0,T  ;  H°  (Rm  ,  B(P))] 


(3'3)  Maj(t>NBi,in     1   K      0    <    t   <    T 


(3.4)      a. (t,x)  is  symmetric 

for    0<t<T,xeR 
—        —  m 
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1         s' 
(3.5)      f  6  L  [0,T  ;  H  (Rm,  P)  ] 

s-1 
n  C[0,T  ;  IT  x (Rm,P) 


(3.6)      uq  e  HS (Rm#P) 


then  (3.1)  has  a  solution  u  such  that 


(3.7)      u6C[0,T;  HS(Rm,P)] 


n  C1[0,T  ;  HS-1(Rm,P)] 


the  solution  is  unique  in  the  larger  class 


C[0,T  ;  H  (Rm,P)]  n  31[0,T  ;  H°(Rm,P)] 


the   following   estimates   hold: 


(3.8,  Hu<t>llr,m£eCKT    (||u0||rrm  +  cPr) 


(3-9)  |||L  u(t)||r.1/m<    c(||f(t)||r.lrm  +   K||u(t)||rrm)     , 


0    <    t    <_    T 

where    0  £   r  £   s' in    (3.8)    and    1    <    r    <    s'  in    (3.9) 
and 

T 
(3.10)  F      =  |  If  (t)  I  I  dt 

o 
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the  constant  c  depends  on  m  and  s. 
Proof.   See  [ 3,  pp.  189  -  194] 
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IV.   General  non-linear  theorem 

The  linear  theorem  is  now  used  to  prove  a  non-linear  existence 
and  uniqueness  theorem  which  is  a  modification  of  [3,  Theorem  II, 
p.  195] .   Instead  of  the  system  Q  [3,  p.  195]  we  consider  the  system 


a      m 
(4.1)       f£  +  Z      G  [u(t)]  iSL_  =  o 
3t    j=1  3  9Xj 


u(0)  =  uq 


(For  present  purposes  the  operator  G  multiplying  r^  can  be  set 

O  d  t 

to  one  and  the  right  hand  side  F  to  zero).   The  quantities  G. (•) 

are  non-linear  operators  that  carry  functions  valued  in 

L1[0,T  ;  HS (R  , p) ]  into  functions  valued  in  L  [0,T  ;  HS (R  ,B(P))] 
m  m 

■  i 

The  difference  here  from  Q   is  that  in  Q   the  operators  G .  act 

pointwise   in  t  (with  t  as  a  parameter)  mapping  functions  on  R 

into  functions  on  R   .   In  (4.1)  the  operators  are  non-local  in 

m 

both  x  and  t  and  act  on  functions  defined  on  [ 0 , T]  xR   .   With 

m 

this  generalization  we  can  apply  the  theorem  to  the  Vlasov-Maxwell 
system  of  plasma  physics. 


Theorem  4.1.   Suppose  s  >  m/2  +  1  .  Let  D  be  a  bounded 

s 
open  subset  of  H  (R  ,P) ,  T  an  arbitrary  positive  number,  and  A(T) 

the  set  L1[0,T  ;  D]  f)  C[  0 ,  T  ;  HS-1  (R  ,P)]  .   Let  G.(-)  be  non-linear 

m  j 

1         s 
operators  on  A(T)  to  L  [  0 , T  ;  H„(R  ,B(P))]  .   Assume  the  following 

*j      m 

conditions  hold: 
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(4.2)       The  quantities  G . [h]  are  bounded  (by  a  constant  A) 
in  the 


H^  -  norm  for  h  6  A(T),  uniformly  in  t  e  [0,T] 


(4.3)       For  each  j  and  t  6  [0,T]  the  map  h(t)  ■+  G.[h(t; 
satisfies  the  condition 


sup  ||G.[h(r)]  -  G.[h(r)]||o^m 


<  y(T)  sup   |h(r)  -  h(r)  |  | 

r<t  °'m 


\i  (T)  -  constant  depending  on  T. 


(4.4)       The  maps  t  ■*■   G.[h(t)]  are  continuous  in  the  H  -norm 


for   each  h  e  A(T) 


(4.5)  For   each   h  e   A(T),    G.[h]  (x)    e    B(P)    is    symmetric    for 

each   t,    xe   [0,T]xR      ,    j    =    1,..,    m 

m 


(4.6)       u   6  D 
o 


Then  there  is  a  unique  solution  of  (4.1)  defined  on  [0,T  ]  ,  where 

i 

0  <  T   <  T  such  that 


(4.7)       u  e  C[0,t'  ;  D]  n  C1[0,t';  HS  1(R  ,P)] 
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Proof.   We  are  guided  as  much  as  possible  by  the  proof  of  that 
for  [3,  Theorem  II,  p.  195]  .   Let  us  assume  for  simplicity  that  D 


is  the  set 


u/u  e  HS(Rm,P)  ,  ||u||     <  R 


and  that   lu       <  R/2  .   The  argument  for  more  general  open  sets 
11  o ' ' s ,m 

D  can  be  done  as  in  [3] .   Let  S  be  the  set  of  functions,  v  ,  on 

i      •  s 

[  0,T  ]  ,  T   <  T  ,  to  H  (R  ,P)  such  that 


v(t)  II     <  R 

1 ' s,m  — 


v(t)  C  C[0,t'  ;  HS-1(Rm,P)] 


The  constant  T   is  to  be  determined  later.   With  the  metric 


v  -  v|  |  |  =  sup   ,|  |v(t)  -  v(t)  |  | 

0<t<T  ' 


S  is  a  complete  metric  space  (see  [3,  Lemmas  2.15,  2.18,  p.  188]) 
If  v  e  S  then  v  6  A(T')  and  we  can  form  the  functions 


(4.8)       g^  =  Gj[v(t)] 


If  we  consider  the  linear  system  (4.1L)  formed  by  replacing  the 
coefficients  in  (4.1)  by  (4.8)  we  see  immediately  that  the  con- 
ditions of  Theorem  (3.1)  are  satisfied,  i.e.  conditions  (3.2),  (3.3), 
(3.4)  are  satisfied  by  virtue  of  (4.4),  (4.2),  (4.5).   Let  w  be 


-32- 


the  solution  to  (4.1,L).   It  follows  by  (3.8),  (3.9)  that  w  has 
the  bounds 


w      <    e       u^  I   _ 
1 ' s,m  —        '  '  o '  ' s,m 


■ 

jt  w     .    <  c  A  e       us,m 
dt   ' ' s-l,m  —         ' '  o ' ' 


For  T   <  In  (2)/cA  then  w  e  S  . 

We  define  the  map  <£  (v)  =  w  from  S  into  itself  and  show  that 
if  T   is  sufficiently  small  <\>    is  a  contraction.   Let 

w  =  <J>  (v) 

w  =  <t>  (v)    ,    v  ,  v  e  s 


Then    (w  -  w)    satisfies   the   system 


a  -  m  a 

f-    (w   -    w)    +      I      G.[v]     §—    (w   -    w) 

at  j=1     D  3x. 


m  ~- 

=    -       E        (G.[  v]     -    G.[  v]  ) 


3=1         J  J  D 


(w   -    w)  (0)    =    0 


Since  v  e  S  then  |  |G_.[  v]  \   ls£/in  1  A  .   since  w  e  S  then  w   6  HS_1  , 
s-1  >  m/2  ,  3 


w      ,    <  R 
1  x . ' ' s-l,m  — 

3 
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It  follows  from  Theorem  (3.1)  and  [3,  Lemma  2.11,  pp.  185  -  186] 
that  for  0  <  t  <  T   <  T 


(w  - 


w) (t)       <  ce 
'    ' 'o,m  — 


cXt 


m 
c   I 
j=l 


G.[  v]  -  G.[  v]  I  I   „   |w   I    ,dt 

y  y       ' 'o, £m' '  x.   s-1 


<  cc  e     R 


m 

E 

j=l 


G.[  v]  -  G.[  v]    -  m    dt 


_  cXT 
<  cce    mRy(T) 


sup 
s<t 


v(s)  - 


v(s^  I  L  m  dt 

1  '  o  ,m 


<  cce    mRy(T)T   sup    |v(s)  -  v(s)|| 

—  i  o ,  m 

s<T 


— •  C  T  — 

<  cce     mRy (T)T     I v  -  v 


c,c  -  Constants 

■ 

Taking  the  supremum  for  t  <_   T   of  the  left  side  of  the  above 

inequality  gives 


1 | | w  -  w | 1 |  <  cceCXT  mRy(T)T'  | | |v  -  v 


For  T   sufficiently  small 


cceC  T  mRy (T)T   <  a  <  1 
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Thus 

|||w  -  w|||  <  a  |||v  -  v||| 

and  4>  is  a  contraction.   Thus  4>  has  a  unique  fixed  point  u   in  S 
which  is  then  a  solution  to  (4.1).   That  u  is  the  unique  solution 
satisfying  (4.7)  follows  from  Theorem  (3.1)  and  the  fact  that  any 
solution  satisfying  (4.7)  is  a  fixed  point  of  <t>  . 
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V.   The  Vlasov-Maxwell  system 

In  this  section  Theorem  4.1  is  used  to  prove  local  existence 
and  uniqueness  of  solutions  to  the  Vlasov-Maxwell  system  (1.1). 
The  assumption  will  be  made  that  the  function  g  has  compact  support. 
By  taking  account  of  the  boundedness  of  the  support  of  the  solution 
it  is  possible  to  put  this  problem  into  a  form  so  that  Theorem  4.1 

applies.   The  proof  gives  solutions  in  a  class  of  Sobolev  spaces 

2 

with  sufficiently  high  order  of  L   differentiability  so  that  by 

Sobolev1 s  lemma  strict  solutions  to  the  problem  are  obtained. 

In  this  section  functions  are  defined  on  R   or  R   and  are 

6      3 

valued  in  R  ,R  ,  or  R   .   Where  it  is  obvious  we  suppress  reference 
to  the  dimension  of  the  range  space  of  functions  in  the  notation, 
i.e.,  HS(R6,R1)  =  HS(R6)  ,  HS(R3,R3)  =  HS(R3)  .   Let  us  assume  for 
the  initial  data  to  (l.la,b)  that 


g  e  HS(R  )  s  >  5 


Eq/Bo  e  HS(R3)  s  >  5 


(5.1) 

supp  g  d  {(x,v)/|x|  <_   R,  |  v  |    R} 


9lls,6  -  M 


E      ,  ,    B      ,  <  N 
o ' ' s, 3    '  '  o  '  ' s, 3  — 


R,M,N  -  constants 


-36- 


Theorem  5.1.   Let  s    5  and  the  initial  conditions  to  (l.la,b) 
satisfy  (5.1).   There  is  a  positive  number  T   and  a  unique  function 


f  6  C[0,T  ;HS(Rc)]n  C1[0,T  ; HS  1(R^)] 
o      b  o        6 


such  that  f  is  a  solution  to  (l.la,b! 


Proof.   Let  the  function  a (v)  be  defined  as  follows 


a(v)  = 


*1iv)     ,    a2(v)  ,  o3(v) 


where 


v±  #  |v|  <  2R 


ai(v)  = 


0  ,  Ivl  >  3R 


a.  (v)  e  C  (R_)  (infinitely  dif f erentiable  function 
with  compact  support) . 


The  system  1.1  is  replaced  with  the  following  system 
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(5.2) 


a)      L.  +  ct(v)-7    f  -    (E  +  a(v)*B)'7,  f  =   0 
o  t  A  v 


f (0)   =  g 


b)     il  -  VxB  =  4?r  ,    ,]  vfdv 

3t  v  J<  2R 


l|+VxE  =  0 


E(0)    =  EQ/    B(0)    =  3Q 


EQ/B0  satisfy  1.1c 


Let   D,A   be    the    following    sets    of    functions 


D   =    {h/h   e    HS(R6,R1)     ,     ||h||  <    2M} 


A    =    L1[0,T;D]     n    c[  0  ,  T;  HS-1  (Rg  ,  R±  )  ] 


Assume  s  j>  5  and  T  an  arbitrary  positive  number. 

For  h  6  A  we  solve  Maxwell's  equation  (5.2b)  with  the  integral 
involving  f  on  the  right  hand  side  replaced  by 
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4ir  vhdv    . 

|  v|<_2R 

Clearly   the    functions 


H^[h]  (x,t)    =    4tt  vihdv      are 


r 

J 

v|<2R 


of    class    L1[  0,T;HS  (R    )]     n   C[  0 ,  T;  HS-1  (R3 )  ]     ,    s    >    5    ,    where    for    each 
t   e   [  0,T] 


(5.3)  l|H.[h]|l       ,    <    C(R)     ||h|l        r         C(R)(2M) 

1  S,J    —  SfO    — 

C (R)  -  a  constant  depending  on  R. 

In  order  to  apply  Theorem  (3.1)  we  write  Maxwell's  equations 
in  the  form  of  a  symmetric  hyperbolic  system. 


Let 


Y    =    [E,B]T    ,    Yo    =    [Eo,Bo]T 


H[h]     =    [  H1,H2,H3,0,0,0]T 


In  terms  of  the  vector  Y  Maxwell's  equations  (5.2b)  with  right 
hand  side  H[h]  takes  the  form 
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3      3- 

(5.4)      §r  Y  +   E   A   g^=  H[h] 
3t      i=1   i  ^xi 


Y(0)  =  YQ 


where  A.  are  symmetric  constant  matrices.   For  h  6  A 


H[h]     6    j}[  0,T;HS(R3,R6)]  f)    C[  0  ,  T;  HS_1  (R3  ,  Rg  )  ] 


where   from    (5.3) 


H[hl  I  I  s,3  -    6C(R)  (2M) 


The   vector   Y      e   HS(R,,R,_)     ,       |Y  ,    <    2N  .      The   constant  matrices 

o  36  ' '    o1 's,3  — 

A.    satisfy 
l 


AiMsJl/3   1    K   =    4^/3)1/2    '    s    >    5 


From  Theorem    (3.1)    we   have    that   the   solution    to    (5.4),    call    it 
Y[  h]     ,    exists    and    is    unique   as    an   element   of 


C[0,T;HS(R3)]    n    C1[0,T;HS"1(R3)] 


and    that 
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Yth]|!s/3(t)  <  ecKt(!!YoMSf3  +c 


f 


H[h]     '       ,  (t)    dt) 


0    <    t    <    T 


It    follows    that 


(5.5) 


Y[  h]  |  |  e    ,  (t)         ecKt  (2N   +    ct6  (C  (R)  )  (2M)  ) 


<    ecKT(2N    +    cT6 (C(R) ) (2M) ) 


and    from    (5.4)    that 


(5.6) 


||Y[h]     -   Y[h]  M0f3(t) 


<    ce 


CKt 


H[h]    -   H[h]  |  I      _  (t)    dt 

o  i  j 


<    6ce         C (R) 


h  -  h||0f6(t)   dt 


t   e    [  0,T]     ,h,h€A 
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Let  E[  h]  (x,t)  ,  B[h]  (x,t)  be  the  fields  obtained  by  solving 
Maxwell's  equations  with  the  right  hand  side  for  formed  with  h  . 
We  form  the  following  non-linear  operators  carrying  functions  in 
A  into  functions  in 


L  [  0,T;H^(Rg,B(R1)  )] 


(elements  of  B(R.)  are  simply  multiplies) 


Gi[  h]  (x,v,t)  =  a±(v) 


Gi+3[  h]  (x,v,t)    =    (E[h]    +   a(v)    x    B[h])i    ,    i   =   1,2,3 


One   now  easily   shows    that   the   operators   G.     ,    — ,    Gfi    satisfy   con- 
ditions    (4.2),    — ,     (4.5).       For   G.     ,    i    =    1,2,3      this    is    obvious 

00 

since  a.  6  C   and  is  independent  of  h  .   For  i=4,-,  6  (4.2) 
is  satisfied  by  virtue  of  (5.5)  taking  into  account  the  bound 
for  ou  (v)  .   Condition  (4.4)  follows  from  the  fact  that 


Y[h]  =  [E[h]  ,  B[h]]T  e  C[0,T;HS(R3)] 


Condition  (4.5)  is  obvious  since  G.  are  scalars.   Condition  (4.3 
follows  from  (5.6)  where 


y  (T)  =  6cc  eCKTC(R) 


c  -  a  constant  depending  on  a (v) 
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It  follows  from  Theorem  4.1  that  the  system  (5.2)  has  a  unique 
solution  f  of  class 

C[0,T  ;D]  n  C1!  0,t'  ;HS_1  (Rg)] 

i 
for  some  0  <  T   <_  T  . 

Let  E[f]  ,  B[f]  be  the  field  associated  with  the  solution  f 
to  system  (5.2).   From  Sobolev's  lemma  (as  stated  for  example  in 

[3,  Lemma  2.1,  p.  182]  one  ascertains  that  E.[f]  ,  B.[  f ]  are  of 

1  ' 

class  C  (R_  x  [0,T  ]   and  are  uniformly  bounded  by  a  constant 

i  i 

A   on  the  interval  [  0,T  ]  .   Using  bounds  on  E  and  B  and  the 

fact  that  the  initial  data  has  support  in  the  set 


{  (x,v)  ,  |x|  <_  R  ,  |v|  <  R}  , 

one  integrates  the  characteristic  equations  for  (5.2)  to  obtain 
bounds  on  the  support  of  the  solution.   The  characteristic  equa- 
tions for  (5.2)  are 


dt 


(5.7) 


g  -   (■  +  iv) 


(as  long  as  |v(t)|    2R  along  the  trajectory).   Here 
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x  = 


,X, |X.|X-) 


V  = 


(v1,v2,v3 


E  = 


(E1,E2,E3) 


B  = 


0    B3  -  B2 


B3   °   "  Bl 


B2  -  Bl   ° 


For  (x,v)  e  R,  let  P(x,v)  =  |v|  .   Let  fi (t)  =  supp  f (t)  -  the  support 

6  R6 

of  f  in  R,  at  time  t.   We  define  the  function 
6 


q(t)  =  sup  P  (x,v) 
fi(t) 


which  bounds  the  support  of  f  in  velocity  space  at  time  t. 
Integrating  (5.7)  for  initial  points  in  the  support  of  the  ini- 
tial data  g  we  have  that 

t 


(5.8) 


v(t)  I  <   v(0)  I  + 


E(t) |  +  |B(t) I |v(t) !  dt 


,t 


<  R  +  6A 


(1  +  q(t))  dt 
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Taking  the  supremum   over  all  such  trajectories  originating  in  the 
support  of  g  gives  the  integral  inequality  for  q(t) 


t 
(5.9)      q(t)  <_    R  +  6>       (1  +  q(t)  dt 

o 

t 
=  R  +  6A  t  +  6A      q(t)  dt 


the  solution  is 


/ .  \  ^    T*    6A  t    ,  6  A  t    ,  > 
q(t)  <_  Re     +  (e     -  1) 


This  gives  a  bound  on  the  support  in  velocity  space  of  the 

i 
solution  f  valid  for  an  interval  of  time  [  0,TQ]  ,  Tq  <_   T   for 

which  q(t)  <    2R  . 

The  time  T   is  obtained  by  solving  the  inequality 

_  6A  t    ,  6A  t    .  >  .  -„ 
Re     +  (e     -  1)  <  2R 


the  solution  is 


t  <  1/6A   In  ((2R  +  1)/(R  +  1))  =  T 


setting 


i     ii 

T   =  min  (T   ,  T  ) 
o 


we  have  that  for  t  e  [  0,T  ] 
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supp  f (t)  C  { (x,v)/jv|    2R} 


Thus  for  t  6  [0,T  ]  the  solution  to  (5.2)  coincides  with  the 

solution  to  (1.1).   It  follows  that  for  the  interval  [  0,T  ] 

o 

the  function  f  is  the  unique  solution  of  class 


C[0,T  ;D]  n  ^[0,?    ;HS_1(R  )]  ,  s  >  5  , 
o  o        6        — 


to  system   (1 .1  )  . 

It  is  an  application  of  Sobolev's  lemma  that  the  solution 
obtained  in  Theorem  5.1  is  a  strict  solution  to  the  problem. 
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